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ABSTRACT 

X-ray studies of normal late-type galaxies have shown that non-nuclear X- 
ray emission is typically dominated by X-ray binaries, and provides a useful 
measure of star formation activity. We have modeled the X-ray evolution of 
late-type galaxies over the ~ 14 Gyr of cosmic history, with an evolutionary 
population synthesis code developed by Hurley et al. Our calculations reveal a 
decrease of the X-ray luminosity-to-mass ratio Ly^/M with time, in agreement 
with observations (Fig. 7a). We show that this decrease is a natural consequence 
of stellar and binary evolution and mass accumulating process in galaxies. The 
X-ray-to-optical luminosity ratio Lx/I/b is found to be fairly constant (around 
~ 10^° ergs^-^Lg Q, Fig. 76), and insensitive to the star formation history in the 
galaxies. The nearly constant value of Lx/L-q is in conflict with the observed 
increase in Lx/I/b from z = Q to 1.4. The discrepancy may be caused by intense 
obscured star formation activity that leads to nonlinear relationship between 
X-ray and B-band emission. 

Subject headings: binaries: close - galaxies: evolution - galaxies: general - stars: 
evolution - X-rays: galaxies - X-ray: binaries - X-rays: stars 
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Introduction 



The X-ray emission of a normal late- type galaxy (i.e., one without an active galactic 
nucle i) is often dominated by the integrated emission of the galacti c X-ray binaries (XRBs) 



(e.g., 



Fabbiano fc White 



2003 



Colbert et a 



be classified into two distinct populations (IBhattacharya fc van den Heuvel 



2004 



Fabbiano 



2006). Galactic XRBs can 



199l|): the 



short-lived (< 10^ yr), high-mass X-ray binaries (HMXBs) and the long-lived (> 10^ yr), 
low-mass X-ray binaries (LMXBs). The X-ray emission from HMXBs is usually regarded to 
trace current star formation because of their short lifetimes, while the X-ray emission from 



LMXBs is more closely re lated to the integrated stellar mass fjPtak et al. 



2003 



2001 



Ranalli et al. 



Grimm et al. 



20031). 



With the observations of ga^ 
Chandra and XMM-Newton (see 



axies at redshift z > 0.1 , either from deep surveys by 



Brandt fc Hasing ei 



analysis of distant galaxy 



2002 



2005 



Nandra et al 



2006 



2002 



ields (2: ^ 0.1 — 4; e.g.. 



Lehmer et al. 



2005 



, for a review) or from stacking 



Brandt et al. 



Georgakakis et al. 



2005 



2003 



200 



Reddy fc Steidel 



Hornschemeier et al. 



2004; 



Laird et al 



20081 ) ■ it has become possible to in yestigate the X-ray 



properties o f normal galaxies at cosmologically significant redshifts flHornschemeier et al 



2000 



2003 



Alexander et al. 



2002 



Georgantopoulos et al. 



2005 



Georgakakis et al. 



2007 



Lehmer et al. 



2006 



2007 



Kim et al. 



2006 



Tzanavaris et al. 



2006 



Rosa-Gonzalez et al. 



20071 ). Previous studies showed that the average X-ray luminosities Lx of normal late- type 
galaxies increase with redshift out to 2; ~ 1.4 — 3, and evolve as (1 + zY'^ ~^ over the redshift 



range ^ ~ — 1.4 for X-ray-detected n ormal galaxies (iNornian et al 



Tzanavaris fc Georgantopoulos 



20081). 



Hornschemeier et al. 



2004 : 



Ptak et al 



2007 



( I2OO2I ) performed a statistical 



X-ray study of spiral galaxies in the Hubble Deep Field-North and its flanking fields using 
the Chandra Deep Field-North 1 Ms data set, and observed a factor of ~ 2 — 3 increase in 
the X-ray-to-optical luminosity ratio Lx/L-q from 2; = to 1.4 for the Le-selected galaxies. 
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b im prove the constraints on the X-ray evolution of late-type galaxies, iLehmer et al 



( 120081 ) studied for the first time how the X-ray properties evolve as a function of optical 



luminosity, stellar mass, and s tar formation rate (SFR) of the galaxies in the Chandra 



Deep Fields North and South (lAlexander et al. 



2003 



Giacconi et al. 



20021 ). It was found 



that there is a significant increase (by a factor of about 5 — 10) in the X-ray-to-optical 
luminosity ratio {Ly^/ L-q) and the X-ray-to-stellar mass ratio (Lx/M) for the galaxy 
populations selected by Lb and M, respectively, over the redshift range of z = — 1.4. 
When analyzing the galaxy samples selected with SFR, these authors found that the 
X-ray luminosity-to-SFR ratio (Lx/SFR) is constant over the entire redshift range for 
galaxies with SFR = 1 — lOOMoyr^^, and that the star formation activity (as traced by 
X-ray luminosity) per unit stellar mass in a given redshift bin increases with decreasing 
stellar mass over the redshift range z = 0.2 — 1, consistent with pr e vious studies on ho w 



star formation act i vity depends on ste . 



Bundy et al. 



2006 



Noeske et al. 



2007a 



la r mass (ICow'ie et al. 



Zheng et al. 



1996 



Juneau et al 



2005 



20071 ). Finally, they extended their 



X-ray analyses to Lyman break galaxies at 2; ~ 3 and estimated that the value of Lx/Lb at 
2; ~ 3 is similar to that at z = 1.4. 

X-ray emission of normal galaxies and its evolution have also been the subject 
of theoretical studies. Using a semi-empirical approa ch to link XRB lifetimes with a 



cosmological evolution of SFR, 



Ghosh fc Whitd (l200l[ ) discussed the imprints left by this 



cosmic SFR on the evolution of X-ray luminosities Lx of normal galaxies. They showed 
that the evolving SFRs can strongly affect the integrated galactic X-ray emission, with 
the possibility of signific ant evolution of the X- ra y luminosities eyen wit hin relatively low 



redshift s z < 1 (see also 



White fc Ghosh 



19981). 



Eracleous et al. 



(120041 ) simulated the 



evolution of X-ray luminosities of XRBs after a burst of star formation with duration of 20 
Myr with a population synthesis method, and found that the 2 — 10 keV luminosity reaches 
a maximum after approximately 20 Myr, and the X-ray luminous phase can be sustained 
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for a period of hundreds of Myr. The results were shown to be insensitive to the initial mass 
function (IMF) and the average mass ratio between accreting and donor stars. However, a 
comprehensive study on the evolution of X-ray populations in galaxies and their relation 
with other properties is still lacking. 

In the present work, we use an evolutionary population synthesis (EPS) code to 
calculate the X-ray luminosity of XRBs and its evolution in a normal late-type galaxy 
over the cosmic history. Meanwhile, we calculate the optical luminosity and the galactic 
mass contributed by stellar populations. The objective of this study is to investigate the 
X-ray evolution of late-type galaxy populations, its dependence on the physical properties 
of galaxies (e.g., optical luminosity, stellar mass and mass-to-light ratio) and on the star 
formation history (SFH), from a theoretical point of view. We will also examine how 
the key parameters, such as IMF, common envelope (CE) efficiency, the binary fraction 
and metallicity, may affect the X-ray emission of the galaxies. In Section 2 we introduce 
the method of calculation and model parameters. In Section 3 the calculated results are 
presented and compared with observations. We summarize in Section 4. Throughout the 
paper, we assume a flat A c old dark matter (AC DM) cosmology with Am q = 0.3, Qa,o = 0.7, 



and Hq = 70 kms ^Mpc ^ (ISpergel et al. 
at z = 1. 



20031 ). which imply a look-back time of 7.7 Gyr 



2. Models 



2.1. Assumptions and input parameters 



We used t he EPS code developed by IHurlev et al 



Liu &: Lil (l2007l . see Appendix A in the paper) and 



(2000 



ZuoetaL 



2OO2I) and updated by 



(I2OO8I ) to calculate the X-ray 



luminosity Lx of XRBs, the optical luminosity Lb and the stellar mass M of the host 
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galaxy, as well as their evolu t ion. T he values of the adopted parameters are the same as the 
default ones in 



Hurley et al. 



( 120021 ) if not mentioned otherwise. 



Previous works (jShapley et al. 



2001 



Persic fc Rephaeh 



2007 



Lehmer et al. 



2008|) have 



already shown that the galactic X-ray emission is closely related to the SF activity. So we 
constructed three cases, i.e., constant SF, star-burst SF and cosmic SF cases, to examine 
their effect. We also examined several key parameters, such as the IMF, CE efficiency 
parameter, the binary fraction and metallicity (listed in Table 1 and discussed below) to 
explore their influence on the X-ray evolution of the galaxies. 

1. Constant star formation case 

In thi s case, we adopt a co nstant SFR = O.25M0yr~^ for stars more massive than 5Mq 



derived by 



Grimm et al. 



(120031 ) in our Galaxy, and take the star formation duration (SFD) 
to be 14 Gyr. For each model, we evolve 10^ primordial binary systems, with the same 
grid of initial parameters (i.e., primary and secondary mass, orbital separation) as Hurley 
et al. (2002). We also evolve 10^ primordial single stars, with initial mass logarithmically 
spaced between O.IM0 and SOM©. In our basic model (i.e., model Ml, listed in Table 1) we 
assume the binary fraction / to be 0.5 and evolve each binary and single star on the grid. 
In the following we describe the assumptions and input parameters in our basic model. 



In order to be in parallel with 



Lehmer et al 



(12001 



(120081 ). we take the IMF of 



Kroupa 



hereafter KROUPAOl) for the mass (Mi) distribution of the primary stars. For the 
secondary stars (of mass M2) and binary orbit, we assume a uniform distribution between 
and 1 for the mass ratio q = M2/M1, and a uniform distribution for the logarithm of 
the orbital separation In a (Hurley et al. 2002). We fix the metallicity to be solar over the 
lifetime of the simulated galaxy. 



We assume that any system entering Roche-lobe overflow becomes circularized and 
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synchronized by tidal interaction between the binary components (Belczynski et al. 2008) 
An important p a rameter in the bin ary evolution is the CE efficiency parameter acE 



(jPaczyhski 



1976 



Iben&Livio 



19931), which describes the efficie ncy of convertin g; orbital 



Zuo fc Li 



2010 



for detail). 



energy into the kinetic energy ejecting the envelope (see §2.1.1 in 
It can often reduce the orbital separation of the surviving binaries by a factor of ~ 100, 
resulting in different outcomes of binary systems. In our basic model, we adopt acE = 0.3, 
which can best model the luminosity function of the galaxy (Zuo et al. 2008). 

We also construct several other models (listed in Table 1) by varying the key input 
parameters described as follows. 

(1) As stated above, variations of the CE parameter can considerably change the 
relative numbers of XRBs. However reliable values of acE are difficult to estimate due 
to lack of understanding of the processes involved, although in the literature it is in the 
range from ~ 0.1 to ~ 3.0 (e.g., Taam & Bodenheimer 1989; Tutukov & Yungelon 1993; 
Podsiadlowski, Rappaport & Han 2003). Here we also adopt acE = 1-0 (Model M2) to 
examine its effect. 

(2) The IMF determines the percentage of high-mass stars, consequently the number 
of XR Bs produced, and the X-ray lum inosity of the galaxy. So we also make use of the 



IMF of 



Kroupa. Tout fc Gilmord ( 11993 



hereafter KTG, model M5), which is much steeper 
in the high-mass end than in KROUPAOl. For the secondary masses (M2), we assume the 
mass ratio q follows a power-law distributi on P(q) oc q", and adopt both the conventio nal 



choice of flat mass spectru m, i.e., a = (IMazeh et al. 



Shatsky fc Tokovinin 



2002 



1992 



Goldberg fc Mazeh 



1994 



our basic model. Ml) and a = 1 (Model M3), since recent data 



are more consistent wit 



( iDalton &: Sarazin 



1 "twins" being a general feature of the close-binary population 



1995 



Kobulnicky &: Fryer 



20071 ) 



(3) Surveys of M dwarfs within 20 pc from the Sun indicate that the binary fraction / 
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may be a function of stellar spectral types ( iFischer &: Marcylll992l ). for examp 



for G stars and f > 0-6 fo r massive 0/B stars in the Cygnus 0B2 association ( iLada 



e, f > 0.5 



Kobulnicky &: Fryer 



20061 : 



20071 ). So we also adopt / = 0.8 (Model M4) for comparison. 

tra-luminous 



(4) Observations of the hosts of XRBs revealed that XRBs, especially n 
X-ray sources (ULXs) may prefer to occur in galaxies with low metallicities ( iMapelli et al. 
20091 ). So we vary the metallicity to examine its effect on the X-ray luminosity evolution by 
taking Z = l.SZ© (Model M6), 0.5Zq (Model M7), O.IZq (Model M8) and 0.02^0 (Model 
M9), in order to compare with our basic model (Ml with Zq). The other parameters in 
these models are the same as the ones in our basic model. 

2. Star-burst star-formation case 

In galaxies like our own Galaxy, continuous star formation processes may last for 
several Gyr, but it is not the case for star-burst galaxies. For example, the Antennae 



f Barnes 


1988; 


Mihos et al. 


1993) 



following combinations of SFD and metallicities: 100 Myr/Z© (Model MIO), 20 Myr/Zo 
(Model Mil) and 100 Myr/O.O2Z0 (Model M12). We assume that the SF is quenched after 
the SFD time and set other parameters to be the same as in our basic model. 

3. Cosmic star-formation case 

With improved observations of star formation processes, cosmic SFH can be constrained 
quite tightly w ithin ~ 30% — 50% up to z ~1 and within a factor of ~ 3 up to 2; ~ 6 



( iHopkins 



2OO4J ). which makes it possible to investigate the cosmic X-ray evolution of 
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galaxies. Here we adopt the derived expression of the SFH in iHopkins fc BeacomI (120061 ) 



PSF[Z) OC < 



1 + 
1 + 

[ (1 + zy'-' 



z < 0.97 

0.97 <z< 4.4J 

4.48 < z, 



and scale the SFR at redshift z = to be the same as that of our Galaxy. Moreover, it is 
known that the cosmic metallicity also evolves strongly with redshift, and galaxies at higher 



redshift tend to have lower metalli c ities (iPettini et al.lll999l: iProchaska et al.ll2003 



20031: 



Kobulnickv &: Kewlev 



2004: 



Kewley fc Kobulnickv 



2007 



Kulkarni et al. 



iRao et al. 



2005 



2007 



Savaglio et al. 



2005 



2009 



Savaglio 



2006 



Wolf fc Podsiadlowski 



So we adopt an empirical 



evolution with 7 = 0.15 ( iKewley &: 



equation Z/Zc:^ oc 10 



-72 



one (KTG93) and a shallower one (IBaldry fc Glazebrook 



Kobulnickv 



2007 



Langer fc Norman 



Peroux et al 



20071 ). 



20061 ) for metallicity 



2007D. We also vary the IMF to a steeper 



2003 



BG03 for short) to examine 



its effect in this case. The other parameters are the same as in our basic model. 



2.2. X-ray luminosity and source type 



We adopt the same proced ure to calculate the 0.5 — 8 keV X-ray luminosities of 
different XRB populations as in lZuo fc Lil ( 120101 ). Mass transfer in XRBs occurs via either 



Roche-l obe overflow or c aptur e of the wind material from the donor star. We use the 



classical 



Bondi fc Hoyld ( 1l944l ) 's formula to calculate the wind accretion rate of the compact 



stars. In the case of Roche-lobe overflow, mass is transferred to the accreting star by way 



of an accretion disk. It is known that accretion dis 
instability if the accretion rate is sufficiently low ( 



<:s in LMXBs are s ubject to the thermal 



van Paradiis 



19961). We discr iminate 



transient and persistent LMXBs according t o the criteria of ivan 



sequence (MS) and red giant donors, and of 



Ivanova &: Kalogeral (|2006[ ) for white dwarf 



ParadiisI ( 119961 ) for main 
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^X,0.5-8keV 



(WD) donors, respectively. The simulated X-ray luminosity is described as follows, 

'7boi'7out-^Edd transients in outbursts, 

r7boimin(Lboi,?7Edd-f'Edd) persistent systems, 
where the bolometric accretion luminosity Lboi — O.lMaccC^ (where Mace is the 
accretion rate and c is the velocity of light), the critical Eddington luminosity 
L-Edd — 47rG'mimpc/(7T = 1-3 x lO^^mi ergs~^ (where ctt is the Thomson cross section, 
rrip the proton mass, G the gravitational constant, and mi the accretor mass in the 
units of solar mass), and r^E d d is the factor t o allow super-Eddington luminosities, taken 



(2) 



to be 5 (jOhsuga et al. 



2002 



Begelman 



2OO2I ) ■ To transform the bolometric luminosity 



i nto the 0.5 — 8 keV X-ra y luminosity, a bolometric correction factor r^boi is introduced 



teelczynski &: Taam 



20041 ). Generally, its value is ~ 0.1 — 0.5 for different types of XRB, 



here we adopt r^boi — 0.1. For transient sources the X-ray luminosity during outbursts 
should be larger than the long-term one by a factor r^out- We take r^out = 0.1 and 1 for the 
short and long-period systems, and the critical periods are adopted to be 1 da y for neutron 
star ( NS) transients and 10 hours for black ho. 



1997 



Garcia et al. 



2003 



Belczynski et al. 



20081 ). 



e (BH) transients, respectively flChen et al. 



2.3. Optical luminosity Lb and stellar mass of the galaxy 

The optical luminosity Lb of a galaxy is mostly from normal stars (both binary and 
single stars). Assume that the stellar radiation can be reasonably approximated as a 
blackbody, the B-band luminosity of a star is calculated with Lb = '^jj^^x ^ "where L is the 
total thermal luminosity of the star, and the radiative intensity Ix = ^hc/\kT^g_i J where 
h is the Planck constant, k the Boltzmann constant, and T^s the effective temperature. 

We also examine the contribution of optical luminosity from accretion disks in XRBs, 
resulting from the reprocessing of X-ray photons. We calculate the optical luminosity Lb 
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from the accretion disk in BH XRBs following iMadhusudhan et al.l ( 120081 ). adopting the 
same temperature profile (i.e., their Eq. [4]) to describe the effective temperature in the 
disk. For NS XRBs we find similar results as BH XRBs. Our calculation reveals that 
optical radiation from accretion disks in XRBs is negligible compared to the overall stellar 
optical luminosity. 



The stellar mass here is the sum of the masses of currently living stars, and does 



not include t 
parallel with 



le contribution frorn compact stars (WDs, NSs, and BHs), in order to be in 



Lehmer et al 



( 120081 ). where they used the rest-frame B — V color and fT-band 



luminosity to estimate the masses of the galaxies. 



3. Results 

3.1. Constant star- format ion case 

Figure 1 shows the calculated values of Lx, Lb, Lx/M, Lx/Lb, Lx/{M/Lb) and M/Lb 
against time in the constant SF case. The five panels from left to right correspond to the 
basic model (Ml), models with acE = 1-0 (M2), a = 1 (M3), / = 0.8 (M4), and KTG93 
IMF (M5), respectively. 

For X-ray luminosities Lx (Fig. la) we plot the contributions from HMXBs and LMXBs 
with dotted and dashed lines, respectively. The X-ray luminosity of HMXBs rises rapidly 
shortly after the first SF, and remains nearly constant afterward, because of continuous, 
constant SF. LMXBs take much longer time to form, and their number is correlated with 
the total star mass, resulting in a long-term increasing trend in X-ray luminosity with 
time. The position of the crossing point of the two lines depends most strongly on the CE 
parameter acE (larger acE results in more LMXBs). 

The optical luminosity Lb (Fig. 16) is contributed by the primary stars (dotted line) 
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and secondary stars (dashed line) in binaries, and by single stars (long-dashed hne). When 
the fraction of binaries is larger (Model M4), more XRBs are produced, leading to larger 
Lx, and hence larger Lx/M, Lx/ and Lx/ {M / L^) ratios compared with those in the 
basic model. The steeper end of IMF (i.e.. Model M5) implies a smaller number of massive 
stars (and smaller Lb), resulting in fewer compact objects and smaller Lx- So the Lx/M 
and Lx/ {M / L-q) ratios in this case are smaller than those in models M1-M4. 

The Lx/M ratio (Fig. Ic) has a clear decreasing trend after the age of tens of Myr. 
Observationally, this phenomena was regarded as the evidence that lower-mass galaxies 
appear to have higher specific star formation rates than massive ones (Brinchmann et al. 
2004; Bauer et al. 2005; Noeske et al. 2007a; Feulner et al. 2005; Zheng et al. 2007). Our 
results show that the decrease of Lx/M is a natural consequence of XRB evolution and 
stellar mass accumulating process in galaxies - the stellar mass always steadily increases 
while the X-ray luminosity changes httle during most of the evolution. 

The Lx/Lb ratio (Fig. Id), similar to Lx, rises in the first several Myr, then remains 
nearly constant afterward. The flattening values are all around 10^° ergs~^L3Q, however 
have severalfold changes among different models. They are caused by the diversity in the 
percentages of both total massive stars and massive stars in binaries, which determine 
Lb and Lx, respectively. The ratios of the peak values of log (Lx/Lb) are roughly 
1 : 2 : 0.5 : 1.5 : 1 from models Ml to M5. 

Figure 2 shows the cumulative X-ray luminosity functions (XLFs) of HMXBs (dotted 
hne) and LMXBs (dashed line) in the top panels. The parameters are the same as in 
Figure 1. Note that HMXBs and LMXBs dominate at relatively high (> 10^^ ergs^^) and 
low (< 10'^^ ergs~^) luminosity ends, respectively. We also show the detailed components of 
XRBs which contribute the total X-ray luminosity separately in the middle (HMXBs) and 
bottom (LMXBs) panels of Figure 2. The dotted, dashed, dash-dotted and dash-dot-dotted 
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lines represent persistent BH XRBs (BHp), transient BH XRBs (BHt), persistent NS XRBs 
(NSp) and transient NS XRBs (NSt), respectively. It is seen that for HMXBs persistent BH 
XRBs contribute most to the XLF; for LMXBs, BH-XRBs (both BHp and BHt) dominate 
the high luminosity end of XLF, while in the low luminosity end transient NS-XRBs play a 
more important role. 

Figure 3 shows the evolution of Lx, Lb, Lx/M, Lx/Lb, Lx/{M/Lb) and M/Lb with 
metaUicities taken to be 1.5Zq, Zq, O.SZ©, O.IZq and 0.02^0, corresponding to models 
M6, Ml, M7, M8 and M9 from left to right, respectively. Note that both Lx and Lb have a 
roughly increasing trend with decreasing metallicity. The values of Lx/M and Lx/ Lb are 
comparable among different models. The values of the peak values are ~ 1 : 1 : 0.7 : 1 : 2 
for log(Lx/M) and ~ 1 : 1 : 0.5 : 1 : 1 for log(Lx/I>B). The corresponding cumulative XLF 
are shown in Fig. 4. 

3.2. Steir-burst star-formation case 

We present the evolution of Lx, Lb, Lx/M, Lx/ Lb, Lx/ {M / Lb) and M/Lb in the 
star-burst case in Fig. 5. Here the metallicity and SFD are assumed to be Zq/IOO Myr (left), 
Z0/2O Myr (middle), and 0.02 Z0/IOO Myr (right). It is seen that the X-ray luminosity 
Lx (Fig. 5a) rises to its peak within the SF episodes (peaked at the age of about 100 
Myr for models MIO and M12, and at about 20 Myr for model Mil), then decreases with 
time, lasting at least several 10^ yr with Lx > 10^^ ergs"^. Note that the small serrations 
emerging in late evolution is mainly caused by LMXBs (i.e.. X-ray transient outbursts due 
to thermal instability in the accretion disks). The optical luminosity Lb (Fig. 56) decreases 
sharply when the SF process stops, since massive stars contribute significantly during the 
star burst episode. The Lx/M (Fig. 5c) (and also Lx/(M/Lb), Fig. 5e) ratio roughly 
follows the trend of Lx- However the slope in this case is much steeper than in Figs. 1 and 
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3, because of the rapid decay of the X-ray luminosity. The L^jL^ (Fig. Sd) ratio in the 
three models are all comparable with those in Figs. 1 and 3, implying that it is intrinsically 
not sensitive to the SFH of the galaxies. 



3.3. Cosmic star-formation case 



Figure 6 shows the evolution of Lx, i^B, Ly^/M, Lx/L- p,, Lx / (M/ Ln) and M/Lu with 



time (left) and redshift z (right), with a cos mic SFH (from 



a cosmic metallicity evolution history (from 



Hopkins &: Beacom 



Langer fc Norman 



20061) and 



20061 ) taken into account. 



Note that the X-ray luminosity Lx is mainly dominated by HMXBs in this case over the 
whole cosmic history because of the enhanced SFR with increasing redshift as a whole. 



This is comparable with the work of 



Lehmer et al. 



(120081), where they found that LMXBs 



on average play a fairly small role in the X-ray emission. The Lx/M ratio has a decreasing 
trend after the age of tens of Myr (or increases with z), similar as in the constant/burst 
SF cases. This also confirms our previous conclusion that the decrease of Lx/M with time 
results from XRB evolution and stellar mass growth in the galaxies. The Lx/ L-^ ratio rises 
rapidly in the first several Myr, then stays around 10^° ergs'^Lgg afterward, similar as in 
the constant/burst SF cases, indicating that the Lx/ L-q ratio is not sensitive to the SFH in 
the galaxies. 

To compare the theoretical predictions with observations, we re-plot Fig. 6 c and Qd 
in Fig. 7a and 75, with redshift z ranging from to 2.0. The solid, dotted and dashed 
lines represent the modeled results with IMFs of Kroupa (2001, KROUPAOl), Kroupa, 
Tout, & Gilmore (1993, KTG93) and Baldry & Glazebrook (2003, BG03), respectively. The 
observational data of log(Lx/M) and log(Lx/i^B) are taken from S hapley et al. (2001, S OI) 



Lehmer et al 



(12008 



for normal late- type galaxies in the local universe (open symbols), 
L08) with the stellar masses 10^°'^ < M/Mq < 10^^'^ (squares. Fig. 7a), B-band luminosities 
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10 



10.5 



< Lb/Lb,q < 1 0^^-^ (filled c i rcles Fig. 76) and 10^°° < Lb/Lb,© < lO^^-^ (filled 



squares, Fig. 76), and 



Zheng et al 



fl2007l . Z07) with 10^° ° < M/Mq < 10^^-^ (diamonds. 
Fig. 7a). In these samples the stellar masses are roughly comparable with our simulated 
ones. To c onvert the SSFR into L-v /M, we have made use of the local Lx— SFR relation 



derived by 



Persic fc Rephaelil ( 120071 ). 



It seems that our simulated log(Lx/M) vs. z relations match the observations quite 
well. Our calculations reveal that the modeled stellar masses are similar to each other 
within 10% when we use different types of IMF, suggesting that the variation of the Lx/M 
ratio is mainly caused by the differences in X-ray luminosity Lx- The values of Lx/M can 
vary by a factor of ~7 between models (KTG93 vs. BG03), as seen in Fig. 7a. 

The nearly constant values of Lx/L^ in Fig. 76 s eem not properly rn atch the observed 



increase in Lx/L-q with z (IHornschemeier et al. 



2002 



Lehmer et al 



20081 ). The discrepancy 



originates from the fact that in our simulations we have roughly Lx oc Lb, giving a fiat 
Lx/L-Q — z relation, while observationally it was found that Lx oc Lg^ (Shapley et al. 
2001; Fabbiano & Shapley 2002; Lehmer et al. 2008), leading to increasing Lx/L-q with 
z. Fabbiano & Shapley (2002) have discussed the X-ray-B-band luminosity correlation, 
and suggested that the nonlinear power law dependency in disk galaxies is likely to be due 
to extinction in dusty star-forming regions, which attenuates light from the B-band more 
effectively than it does in the X-ray band. This means that the intrinsic B-band luminosity 
calculated here is generally larger than the measured one, which suffers local extinction in 
the galaxies. Supporting evidence for this hypothesis also comes from the strong positive 
correlations between Lx/L-q and the ultraviolet d ust-extinction meas ure (Lir + Luv)/Luv, 
and the correlation between Lx/ L-q and IR color (ILehmer et al.ll2008[ ). Thus, if the increase 
in Lx I L -R with z in the late- type galaxies is due to an increase in their star formation 



activity (IFabbiano &: Shapley 



2002 



Lehmer et al. 



20081 ). our results are compatible with 
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the observational data at least qualitatively. 

Figure 8 shows the corresponding cumulative X-ray luminosity functions. In general 
they are similar to those in the case of constant SF. 

Our simulations are obviously subject to many uncertainties. For current population 
synthesis investigations it is difficult to tell confidently which parameter combinations 
are the best or most realistic by comparison with observations, since there are many 
uncertainties in the (both explicit and implicit) assumptions and input parameters, and 
simplifications in the treatment of the detail evolutionary processes. For example, the 
simulated X-ray luminosity Lx depends on the adopted values of several parameters, such 
as the bolometric correction factor ?7boi, the common envelope efficiency parameter and 
so on, which may alter its value severalfold. This further affects the values of Lx/M and 
Lx/Lb since the stellar mass M and the optical luminosity Lb is not sensitive to these 
parameters. Despite these limitations it has become possible to investigate the overall 
evolution of stars in galaxies with population synthesis, and to draw useful information by 
comparing theoretical predictions with observations. 



4. SUMMARY 



We have used an EPS code to calculate the X-ray evolution of late-type galaxies, to 
investigate the relations between the X-ray luminosity and other physical properties (i.e., 
optical luminosity, stellar mass, and star formation history) of the galaxies, and how these 
relations are influenced by the input parameters of star formation and evolution (e.g., 
SFH, IMF, metallicity, and common envelop efficiency, etc). The res ults are compared 



with multi-wavelength analyses and o bservations of late-type galaxies (iShapley et al 



Zhene et al. 


2007; 


Lehmer et al. 


2008) 



2001 



20081 ). In different cases of SF, we find a common feature of 
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decreasing X-ray luminosity-to-mass ratio L^/M with time, in agreement with observations 
(Fig. 7a). We show that the decrease of Lx/M results from slow evolution of Lx of XRBs 
and the stellar mass accumulation with time in galaxies, without requiring that lower mass 
galaxies have higher SSFR than more massive ones suggested before (Brinchmann et al. 
2004; Bauer et al. 2005; Noeske et al. 2007a; Feulner et al. 2005; Zheng et al. 2007). The 
Lx/ Lb ratios in all cases rise rapidly in the first ~ 10® yr to ~ 10'^° erg s~^Lgg, then 
stay nearly constant afterward for a given model, and are not sensitive to the SFH details 
in the galaxies (see Fig. 7b). This seems t o be in conflict with the observed increase in 



2002 



Fabbiano &: Shapley 



2002 



Lehmer et al. 



20081) 



Lx/ L^ with z (iHornschemeier et al 
The discrepancy may be due to different obscured star formati on activities in galaxies at 



higher redshifts (IFabbiano fc Shapley 



2002 



Lehmer et al. 



20081 ). This will be investigated 



by future high-resolution X-ray and optical observations of galaxies at high redshifts. 
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Table 1: Parameters adopted for each model. Here acE is the CE efficiency parameter, q 
the initial mass ratio, IMF the initial mass function, / binary fraction, Z metallicity in solar 
units, an d SFD the duration of star formation in the simulated galaxy. 





Model 




P(q) 


IMF 


/ 


Z 


SFD 














Ze 


(Myr) 




Ml 


0.3 


oc 


KROUPAOl 


0.5 


1.0 


14000 




M2 


1.0 


oc g° 


KROUPAOl 


0.5 


1.0 


14000 




M3 


0.3 


oc q^ 


KROUPAOl 


0.5 


1.0 


14000 




M4 


0.3 


cx g° 


KROUPAOl 


0.8 


1.0 


14000 


Constant SF 


M5 


0.3 


oc g° 


KTG93 


0.5 


1.0 


14000 




M6 


0.3 


ocgO 


KROUPAOl 


0.5 


1.5 


14000 




M7 


0.3 


ocgO 


KROUPAOl 


0.5 


0.5 


14000 




M8 


0.3 


oc g° 


KROUPAOl 


0.5 


0.1 


14000 




M9 


0.3 


oc g° 


KROUPAOl 


0.5 


0.02 


14000 




MIO 


0.3 


oc g° 


KROUPAOl 


0.5 


1.0 


100 


Star-burst SF 


Mil 


0.3 


oc g° 


KROUPAOl 


0.5 


1.0 


20 




M12 


0.3 


oc g° 


KROUPAOl 


0.5 


0.02 


100 
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Fig. 1. — The evolution of (a) the X-ray luminosity (Lx), (b) the optical luminosity (L^), 
(c) the X-ray luminosity-to-stellar mass ratio (Lx/M), (d) the X-ray-to-B band luminosity 
ratio (Lx/L-q), (e) the Lx/{M/L^) ratio, and (/) the stellar-mass-to-B band luminosity 
ratio {M/Lb) with time in the constant star formation case. Here the metallicity is Zq. We 
assume that the secondary mass distribution follows the power-law P{q) = q°' and the binary 
fraction is /. The left panels show the results in the basic model with acE = 0.3, a = 0, 
/ = 0.5 and KROUPAOl IMF. The other models from left to right are with acE = 1-0, 
a = 1, / = 0.8 and KTG93 IMF, respectively. 
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Fig. 2. — The cumulative X-ray luminosity functions (Top: HMXBs + LMXBs; Middle: 
HMXBs; Bottom: LMXBs). The model parameters are the same as in Fig. 1. The dotted, 
dashed, dash-dotted and dash-dot-dotted lines in middle and bottom panels represent black 
hole persistent (BHp) and transient (BHt), neutron star persistent (NSp) and transient (NSt) 
sources, respectively. 
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Fig. 3. — Same as Fig. 1 but for different metallicities, which are taken to be 1.5Zq, Zq, 
0.5^0, O.IZq, and 0.02Zq from left to right, corresponding to models M6, Ml, M7-M9, 
respectively. 
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Fig. 4. — The cumulative X-ray luminosity functions (Top: HMXBs + LMXBs; Middle: 
HMXBs; Bottom: LMXBs). The model parameters are the same as in Fig. 3. 
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Log(oge/Myr) 

Fig. 5. — The Lx, Lb, Lx/M, Lx/Lb, Lx/{M/Lb), M/Lb evolution with time in the star- 
burst case. Here the metalhcities and SFH are Zq/100 Myr (left, MIO), Zq/20 Myr (middle, 
Mil), and O.O2Z0/IOO Myr (right, M12), respectively. 
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Fig. 6.— The Lx, Lb, Lx/M, Lx/Lb, Lx/{M/Lb), M/Lb evolution with time (left) and 
redshift z (right), respect ively. Here we have assumed a cosmic star form ation history (from 



Hopkins fc Beacomll2006l ) and a cosmic metallicity evolution history (from 

2mm- 



Langer fc Norman 
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Fig. 7. — Same as panels c and c? in Fig. 6 but enlarged with z = — 2.0 for comparison 
with observations. The solid, dotted and dashed line represent modeled results with IMFs 
of Kroupa (2001, KROUPAOl), Kroupa, Tout & Gilmore (1993, KTG93), and Baldry & 
Glazebrook (2003, BG03), respectively. Also shown are the measured values of log(Lx/M) 
(left panel) selected by M (lO^^-^ < M/Mq < 10^^ ^ squares), and log(Lx/LB) (right panel) 
selected by Lb (lO^^-^ < Lb/Lb,q < lO^^-\ filled circles; 10^°-° < Lb/Lb,© <iO^^:i^filled 



squares), respectively, for stacked normal late- type galaxy samples derived by 



Lehmer et al. 



(12008 



, LO S). The syrnbols a re the same as in Fig. 10 in LOS. We converted the SSFR to 



Lx/M in 



Zheng et al. 



( l2007l . Z07) samples of the corresponding stellar mass bin (10 



10.0 



< 



M/Mq < 10^"-^ diamonds), using the local Lx— SFR relation derived by lPersic &: Rephaelil 



(120071 ). The data for normal late-type galaxies in the local universe (open symbols) are from 
the Shapley et al. (2001, SOI) samples. 




Fig. 8. — The cumulative X-ray luminosity function (Left: H-|-L; Middle:H; Right: L ). Here 
we have assumed a cosmic star formatio n history (from Hopkins fc Beacom 20061 ) and a 
cosmic metallicity evolution history (from lLanger fc Normanll2006l ). 



